Microscopic Formulation of Interacting Boson Model for Rotational Nuclei 
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We propose a novel formulation of the Interacting Boson Model (IBM) for rotational nuclei with 
axially-symmetric strong deformation. The intrinsic structure represented by the potential energy 
surface (PES) of a given multi-nucleon system has a certain similarity to that of the corresponding 
multi-boson system. Based on this feature, one can derive an appropriate boson Hamiltonian as 
already reported. This prescription, however, has a major difficulty in rotational spectra of strongly 
deformed nuclei: the bosonic moment of inertia is significantly smaller than the corresponding 
nucleonic one. We present that this difficulty originates in the difference between the rotational 
response of a nucleon system and that of the corresponding boson system, and could arise even if 
the PESs of the two systems were identical. We further suggest that the problem can be cured by 
implementing L-L term into the IBM Hamiltonian, with coupling constant derived from the cranking 
approach of Skyrme mean-field models. The validity of the method is confirmed for rare-earth and 
actinoid nuclei, as their experimental rotational yrast bands are reproduced nicely. 

PACS numbers: 21.10.Re,21.60.Ev,21.60.Fw,21.60.Jz 



The atomic nucleus is a strongly interacting many- 
body quantal system which has collective properties re- 
sulting in various deformed shapes. If a nucleus is 
strongly deformed, it rotates, exhibiting characteristic ro- 
tational band structure with remarkable regularity. Such 
rotational motion can be viewed as a manifestation of 
symmetry restoration mechanism of Nambu, and has at- 
tracted much attention in nuclear physics from various 
viewpoints [J, Q. 

The Interacting Boson Model (IBM) 0, [| has been 
successful in phenomenological studies for describing low- 
lyingquadrupole collective states of medium-heavy nu- 
clei [3]. The major assumption of the IBM is to employ 
L = + (s) and L — 2 + (d) bosons which reflect the 
collective S and D pairs of valence nucleons [|| . The mi- 
croscopic foundation of the IBM has been studied exten- 
sively so as to derive an IBM Hamiltonian starting from 
nucleon degrees of freedom ■ A new approach to de- 
rive the IBM Hamiltonian has been presented recently 
Q. In this approach, the potential energy surface (PES) 
with quadrupole degrees of freedom, obtained from the 
mean-field calculation with the Skyrme Energy Density 
Functional (EDF) 0, [l(j , is compared to the correspond- 
ing PES of IBM to obtain the parameters of IBM Hamil- 
tonian. As a Skyrme EDF gives universal descriptions of 
various nuclear properties [1, l9Hl~l| , one can derive the 
IBM Hamiltonian basically for all situations in a unified 
way. This method turned out to be valid particularly for 
nuclei with weak to moderate quadrupole deformation, 
and has been practiced extensively [l2[ • When a nucleus 
is well deformed, however, the nucleonic rotational spec- 
trum appears notably and systematically different from 
the corresponding bosonic one, being manifested by too 
small bosonic moment of inertia as compared to the cor- 
responding fermionic one [1, [l^l ■ 

This kind of difference has been known as a result of 



limited degrees of freedom, s and d bosons only, in many 
cases [H, UH ■ In order to remedy this problem, another 
type of nucleon pairs, e.g., L = 4 + (G) pair, and the cor- 
responding boson image (g boson) have been introduced, 
and their effects were renormalized into sd boson sector, 
yielding IBM Hamiltonians consistent with phenomeno- 
logical ones [til. ItI [l3l - l2lj . In the mean time, the validity 
of IBM for rotational nuclei was analyzed in terms of the 
Nilsson model [22|, coming up with the criticism that 
the iS-D-pair truncation may be far from sufficient for de- 
scribing intrinsic states of strongly deformed nuclei, and 
this naturally casts a question on the applicability of the 
IBM to rotational nuclei in particular. While it has been 
reported that the SD-pa.ii dominance holds to a good 
extent in intrinsic states of rotational nuclei [HI, 24 1 , 
there has been no conclusive mapping procedure from nu- 
cleonic systems to IBM ones covering rotational nuclei. It 
is thus of much interest to revisit this issue with the newly 
proposed method of Ref . @ , looking for a prescription to 
cure the afore-mentioned problem of too small moment 
of inertia. 

In the method of Ref. Q , we calculate the energies of 
nucleonic and bosonic intrinsic states representing vari- 
ous shapes, and obtain PESs. We then determine param- 
eters of the IBM Hamiltonian so that the bosonic PES 
becomes similar to the nucleonic one [|[. These intrin- 
sic states are at rest with rotational frequency uj—0. In 
this paper, we move on by one step further with non-zero 
rotational frequency w ^ 0. Actually we analyze the re- 
sponses of the nucleonic and bosonic intrinsic states by 
rotational cranking with infinitesimal oj. From such re- 
sponses, one can extract the most important rotational 
correction to the IBM Hamiltonian. 

The nucleon intrinsic state \4>f) is obtained from the 
Hartree-Fock plus BCS (HF+BCS) calculation. Skyrme 
SkM* interaction [26[ is used throughout, while different 
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FIG. 1: (Color online) Overlap between the intrinsic state and 
its rotation at angle fi for 146_154 Sm and 23 °- 238 TJ nuclei for 
(a)(c) fermion (HF+BCS) and (b)(d) boson (IBM) systems. 



Skyrme forces do not alter the conclusion. 

For boson system, we consider the IBM-2, because it is 
closer to a microscopic picture than the simpler version of 
IBM. The IBM-2 is comprised of proton L = 0+ (s n ) and 
2+ (dn) bosons, and of neutron L — + (s v ) and 2 + (d v ) 
bosons 0. We take the standard IBM-2 Hamiltonian, 



Hb = en d + KQ n ■ Q v , 



(1) 



where rid — + n<iv with fi^ (rid-*) being the proton 
(neutron) d-boson number operator and Q p = s' p d p + 

d p Sp + Xp[dpdpY 2 ^ ■ Here e, k, and X-n,v are parameters, 
and their values are determined by comparing nucleonic 
and bosonic PESs following Refs. [1, [llj]. The boson 
intrinsic state \4>b) is written in general as a coherent 
state (H EH 
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where |0) and a pp represent the boson vacuum (inert 
core) and amplitude, respectively. 

We now look into the problem of rotational response. 
We shall restrict ourselves to nuclei with axially symmet- 
ric strong deformation, because this problem is crucial to 
those nuclei but is not so relevant to the others. An axi- 
ally symmetric intrinsic state is invariant with respect to 
the rotation around the symmetry (z) axis. This means 
a Pti = for fi ^ in cq. @ in the case of bosons. Such 
intrinsic states of nucleons and bosons are supposed to be 
obtained as the minima of the PESs. Let us now rotate 
the axially symmetric intrinsic states about the y-axis by 
angle fi. Figure Q] shows the overlap between the intrinsic 



state \4>x) and the rotated one \4>' x ) = e~ 
X stands for either fermion (X—F) or 
Here L y denote the y-component of the angular momen- 
tum operator. We take 146 ™ 154 Sm and 23 °- 238 u as ex- 
amples. Some of these nuclei are good examples of SU(3) 
limit of IBM [25|. 

Figures [Ua,c) an d Hb,d) show the overlaps for nu- 
cleons and bosons, respectively. For Sm isotopes, the 
parameters of Hb are taken from [l2j], while the param- 
eters for U isotopes are determined as e w 0.100 MeV, 
k w —0.18 MeV, and Xw w Xv w —1-0. These parameters 
are used throughout this paper. 

In each case, the overlap is peaked at fi — 0° with 
the value unity, and decreases with fi. The nucleonic 
overlaps are peaked more sharply, whereas boson ones 
are damped more slowly. It is clear that as a function of 
fi, boson rotated intrinsic state changes more slowly than 
the corresponding nucleon one, due to limited degrees of 
freedom for bosons. 

We point out that the overlap becomes narrower in 
fi with the neutron number N for Sm isotopes (see 
Fig. [Ua,b)). This is related to the growth of deforma- 
tion. On the other hand, there is no notable change in the 
overlap for these U isotopes, because pronounced prolate 
minimum appears always at fii ~ 0.25 in their PES. 

The nucleon-boson difference of the rotational response 
discussed so far suggests that the rotational spectrum of 
a nucleonic system may not be fully reproduced by the 
boson system determined by the mapping method of Ref. 
8] using the PESs at rest. In fact, it will be shown later 
that the moment of inertia of a nucleon system differs 
from the one calculated by the mapped boson Hamilto- 
nian. We then introduce a term into the boson Hamil- 
tonian, so as to keep the PES-based mapping procedure 
but incorporate the different rotational responses. This 
term takes the form of L ■ L where L denotes the boson 
angular momentum operator. This term is nothing but 
the squared magnitude of the angular momentum with 
the eigenvalue L(L + 1), and changes the moment of in- 
ertia of rotational band keeping their wave functions. A 
phenomenological term of this form was used in the fit- 
ting calculation of IBM, particularly in its SU(3) limit 
25] , without knowing its origin or physical significance. 

We adopt, hereafter, a Hamiltonian, H' B , which in- 
cludes this term with coupling constant a : 

H' B = H B +aL- L, 



(3) 



where Hb is given in Eq. ([T]) . The aL ■ L term will be re- 
ferred to as LL term hereafter. The LL term contributes 
to the PES in the same way as a change of d-boson en- 
ergy Ae = 6a (see eq. (QJ), because the PES at rest (i.e., 
lu=0) is formed by the boson intrinsic state \4>b) contain- 
ing no d±\ component. Hence, by shifting e slightly, we 
obtain the same PES as the one without the LL term, 
and consequently the other parameters of mapped Hb 
remain unchanged. 

We now turn to the determination of a in Eq. ([3]). 
First, we perform the cranking model calculation for the 
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fermion system to obtain its moment of inertia, denoted 
by Jp, in the usual way B. By taking the Inglis-Belyaev 
(IB) formula, we obtain [29|, |3Cj 
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where energy Ei and v-factor Vi of quasi-particle state 
% are calculated by the HF+BCS method of Ref. [3lj |. 
Here, is the nucleon angular momentum operator, and 
k means the axis of the cranking rotation, being either x 
or y, as z-axis. Based on the earlier argument, the y-axis 
is chosen. 

Next, the bosonic moment of inertia, denoted as Jb, is 
calculated by the cranking formula of Ref. [HJ with d±i 
being mixed, to an infinitesimal order, into the coherent 
state 1 4> B ) in Eq. ([2j: 



Jb = lim — 
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where U) is the cranking frequency, a p ±\ denotes the am- 
plitude for d p ±\, and stands for the boson angular 
momentum operator. Note that a p ±i oc U) at this limit, 
leading Jb to a finite value. 

The value of a is determined for individual nucleus so 
that the corresponding bosonic moment of inertia, Jb in 
eq. ([5]) becomes equal to J F in eq. This prescrip- 
tion makes sense, if the nucleus is strongly deformed and 
the fixed intrinsic state is so stable as to produce indi- 
vidual levels of a rotational band through the angular 
momentum projection in a good approximation. The re- 
sultant excitation energies should follow the rotor for- 
mula E x oc L(L + 1) for L being the angular momentum 
of the level. The present prescription with the LL term 
should be applied only to certain nuclei which belong to 
this type. We introduce a criterion to select such nu- 
clei in terms of the ratio -R4/2 = E x (Af)/E x (2f), and 
set a minimum value for this. Empirical systematics [33[ 
suggests that the evolution towards stronger deformation 
continues as the number of valence nucleons increases, 
but this evolution becomes saturated beyond -R4/2 ~ 3.2. 
Namely, for the nuclei with -R4/2 > 3.2, the deformation 
is considered to be evolved sufficiently well, and we take 
-R4/2 > 3.2 as the criterion to apply the LL term. This 
discrete criterion is also for the sake of simplicity, but 
the major discussions of this work do not depend on its 
details. 

Figures [2a) ~ (c) show the moments of inertia for Sm 
and U isotopes. In these figures, Jb calculated with the 
LL term (w/ LL), Jb calculated without it (w/o LL), 
and Jf are compared. Experimental ones determined 
from the 2+ levels j3f| are shown also. 

We divide Sm isotopes into two categories according 
to the criterion defined above. First, the ratio -R4/2 is 
calculated without the LL term, leading to 152 ~ 158 Sm 
with i? 4/2 > 3.2 and 146 - 150 Sm with i? 4/2 < 3.2. For the 
former category, the LL term should be included, and 
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FIG. 2: (Color online) (Upper panels) Moments of inertia in 
the intrinsic state for (a) 146 - 150 Sm, (b) 152 - 158 Sm and (c) 
232-236 Ui ca i cu i atec i b y IBM with (w/) and without (w/o) 
the LL term and by IB formula. Experimental data taken 
from _Ex(2^~) [35[ are also shown. (Lower panel) The derived 
a value for (d) 146 - 150 Sm, (e) 152 " 158 Sm and (f) asa-ase^ 



Fig. [2Kb) demonstrates that the LL term produces sig- 
nificant effects so as to be just enough for the agreement 
to experiment. To be more precise, the experimental 
value becomes large for N = 90, and looks nearly flat 
for N > 92, being 35 ~ 40 MeV^ 1 . The enlargement 
of moment inertia means that the value of a is negative. 
The IB formula reproduces this trend quite well, which 
is inherited to bosons by the present method. 

Although the LL term should not be used for the cate- 
gory depicted in Fig. [^a), we shall study some features. 
The increase of Jf and Jb with increasing N can be 
seen. Although experimental moment of inertia exhibits 
a gap between Figs. [2^ a) and (b) (from N —88 to 90), 
neither Jb nor Jp does not follow this trend, showing 
only gradual changes. This could be due to the absence 
of the particle number conservation in the Skyrme EDF 
calculation. We do not touch on this point in this paper. 

Figures EJd) and (e) show, respectively, the derived 
a value for 146_150 Sm and 152_158 Sm. First we notice 
an overall trend that a does not change so much, while 
the IB value of Jf changes by an order of magnitude. 
Although the a values for 146_150 Sm do not make much 
sense, this is of certain interest. 

FigurelDJc) shows the moments of inertia for 232 - 236 tj j 
which are rather flat. We point out that the calculated 
moment of inertia, J f =Jb with the LL term, turned 
out to be about twice large as that of 152_158 Sm. This 
dramatic change is consistent with experiment, although 
somewhat overshoots experimental changes. 

We shall discuss eigenvalues of H' B in eq;_j3]) obtained 
by the diagonalization by NPBOS code jllfT We first 
investigate to what extent E x (2^) is lowered by the LL 
term. Figure [3] shows the fraction of this lowering, by 
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FIG. 3: (Color online) Fraction of _E X (2+) with the LL term, 
normalized with respect to E x (2f) without the LL term, for 
(a) Sm and (b) U isotopes. Arrows represent the LL matrix 
elements normalized by i? x (2+) (w/o LL). 
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FIG. 4: (color online) Experimental [3g] (symbols) and calcu- 
lated (curves) yrast spectra for (a) Sm and (b) Gd isotopes. 



normalizing it with respect to the E x (2^) without the 
LL term, for (a) Sm and (b) U isotopes. This lowering 
is, as indicated by the arrows in Fig.[3j >30% for 92 ~ 96 Sm 
and >60% for 140 ~ 144 u. On the other side, it is almost 
vanished or quite small for N = 84 ~ 90. Thus, it may 
not affect the IBM description much, even if one keeps 
the LL term in all nuclei. We do not take it, because the 
present derivation does not give physical basis for the LL 
term for nuclei without strong deformation. 

Figure @] shows the evolution of low- lying yrast spectra 
for (a) Sm and the neighboring (b) Gd isotopes as func- 
tions of N. For both Sm and Gd isotopes, the LL term 
is included for N ^ 90, but is not included for N ^ 88, 
based on the criterion discussed above. The IBM param- 
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FIG. 5: (Color online) Level schemes of (a) 154 Sm, (b) 156 Gd, 
(c) 230 Th and (d) 232 U nuclei. Calculated spectra with (w/) 
and without (w/o) the LL term, are compared with experi- 
mental [33| spectra. 



eters for Gd isotopes are derived similarly to those used 
for Sm isotopes. Figures. |4ja) and EJb) indicate that 
calculated spectra become more compressed with TV and 
exhibit rotational feature for N ^ 90, similarly to the ex- 
perimental trends [35[ . One notices a certain deviation at 
N = 88, where the Skyrme PES favors stronger deforma- 
tion and the calculated excitation energies are somewhat 
too low (T^ |. 

Figure[S] shows yrast levels of 154 Sm, 156 Gd, 230 Th and 
230 U nuclei as representatives of rotational nuclei. The 
LL term is included for these nuclei, as they fulfill the 
criterion. For 230 Th, the parameters of Hb take almost 
the same values as those for the 232 U nucleus. A nice 
overall agreement arises between the theoretical and the 
experimental 35] spectra, and the contribution of the 
LL term to it is remarkable. Particularly for 154 Sm and 
230 Th, the calculated spectra look nearly identical to the 
experimental ones. 

We comment on side-band levels. The deviations of 
/3-bandhead (0^) and 7-bandhead (2 J) energies are im- 
proved by tens of keV by the LL term. However, these 
band- head energies are still much higher than experimen- 
tal ones. Thus, there are still open questions on side-band 
levels. On the other hand, the relative spacing inside the 
bands is reduced by hundreds of keV, producing certain 
improvements. 

We mention some studies deriving a collective Hamilto- 
nian from a given EDF where the mean-field PES supple- 
mented with zero point rotational and vibrational correc- 
tions are treated as an effective potential [36- 38] . A gen- 
eralized kinetic energy term emerges in such approaches. 
In the present work, we compare the results of Skyrme 
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EDF with the corresponding results of the mapped boson 
system, at the levels of the PES and the rotational re- 
sponse. The kinetic energies of nucleons are included in 
both levels, while the rotational kinetic-like boson term 
appears from the latter. 

In summary, we have proposed a novel formulation of 
the IBM for rotational nuclei. The rotation of strongly 
deformed multi-nucleon system differs, in its response to 
the rotational cranking, from its boson image obtained 
by the mapping method of Ref. [|[ where the PES at 
rest is considered. Significant differences then appear 
in moment of inertia between nucleon and boson sys- 
tems. We have shown that this problem is remedied 
by introducing the LL term into the IBM Hamiltonian. 
The effect of the LL term makes essential contribution 
to rotational spectra, solving the longstanding problem 
of too small moment of inertia microscopically. Experi- 
mental data are reproduced quite well, without any phe- 
nomenological adjustment. The mapping of Ref. [8j ap- 



pears quite sufficient for vibrational and 7-unstable nu- 
clei, and the present study makes the IBM description of 
strongly deformed nuclei sensible theoretically and em- 
pirically. Thus, we seem to have come to the stage of 
having microscopic basis of the IBM in all situation at 
the lowest order. On the other hand, this achievement is 
partly due to the successful description of Skyrme mean- 
field model. The feature discussed in this paper is related 
to the question as to whether the IBM can be applied to 
deformed nuclei or not [22| . The present work indicates 
that the rotational response is substantially different be- 
tween fermions and bosons, but the difference can be 
incorporated into the IBM in a microscopic way. 
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